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We study a non-Hermitian version of the Rabi model, where a two-level system is periodically 
driven with an imaginary-valued drive strength, leading to alternating gain and loss. In the Floquet 
picture, the model exhibits PP symmetry, which can be broken when the drive is sufficiently strong. 

We derive the boundaries of the PP phase diagram for the different resonances by doing perturbation 
theory beyond the rotating-wave approximation. For the main resonance, we show that the non- 
Hermitian analog of the Bloch-Siegert shift corresponds to maximal PP breaking. For the higher- 
order resonances, we capture the boundaries to lowest order. We also solve the regime of high 
frequency by mapping to the Wannier-Stark ladder. Our model can be experimentally realized in 
waveguides with spatially-modulated loss or in atoms with time-modulated spontaneous decay. 


I. INTRODUCTION 

In recent years, there has been substantial interest in 
systems that are symmetric under parity-time reversal 
{PT). A non-Hermitian Hamiltonian can still have real 
eigenvalues if it is VP symmetric [1]. In general, as a pa¬ 
rameter is varied, a T^T-symmetric Hamiltonian under¬ 
goes a transition from a VP -symmetric phase (real eigen¬ 
values) to a T^T-broken phase (complex eigenvalues). In 
the T^T-symmetric phase, the system exhibits periodic 
oscillations in time, while in the T^T-broken phase, it 
exhibits exponential growth. This discovery opened up 
the possibility of studying what new physics arises in 
non-Hermitian systems |2HT6]. 

While research on VP symmetry has focused on time- 
independent Hamiltonians, several recent works have 
studied time-dependent Hamiltonians [IMS]. A recent 
work considered a non-Hermitian Rabi model, where a 
two-level system is periodically driven in time with an 
imaginary-valued drive strength, i.e., there is alternating 
gain and loss [18]. There is a rich VP phase diagram as 
a function of drive frequency and strength. Furthermore, 
an imaginary drive can induce Rabi oscillations between 
the two states [IS] . 

In this paper, we use perturbation theory to derive 
the boundaries of the VP phase diagram for the non- 
Hermitian Rabi model. In particular, we capture higher- 
order effects beyond the rotating-wave approximation. 
The model is defined as 
UJo 

H = + ‘iiXcFx coswt, (1) 

which represents a two-level system with energy differ¬ 
ence Wo, driven periodically at frequency w. The drive 
strength i\ is imaginary, leading to periodic gain and loss 
in time. (Jx and cr^ are Pauli operators. The goal is to 
calculate, as a function of parameters, whether evolution 
under Eq. 0 leads to periodic dynamics (PT-symmetric 
phase) or to exponential growth (T^T-broken phase). In 
the Floquet picture [20j|2l], we want to see whether the 
Floquet quasienergies of Eq. 0 are real or complex. The 
strategy here is to recall the known results for the Her- 



FIG. 1. (Color online) PP phase diagram, showing 
"PT-symmetric (white) and "PT-broken (gray) phases. Also 
shown are the boundaries for the single-photon resonance (red 
dashed line) and Bloch-Siegert shift (blue dotted line) calcu¬ 
lated using perturbation theory at next-to-leading order. 

mitian Rabi model (with a real-valued drive strength), 
and then let the drive strength become imaginary and 
see when the Floquet quasienergies become complex. 

We first derive the boundary for the “single-photon” 
resonance (Fig. [^. In particular, we find the curva¬ 
ture due to counter-rotating terms and identify the non- 
Hermitian analog of the Bloch-Siegert shift [^. Then we 
derive the boundaries of the “multi-photon” resonances 
to lowest order (Fig.[^ and show that the width of succes¬ 
sive resonances decreases exponentially. We also derive 
the boundaries in the limits of low and high frequencies, 
where for the latter, we use a modified high-frequency 
Floquet analysis. 

Equation ([^ can be implemented in experiments with 
waveguides, where propagation in space corresponds to 
propagation in time [23]. Recent waveguide experiments 
have implemented the Hermitian Rabi model |24j . By 
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FIG. 2. (Color online) Zoomed-in view of VT phase diagram, 
showing multi-photon resonances at cj « uJol{‘2n -\- 1). Also 
shown are the lowest-order boundaries for some multi-photon 
resonances (red dashed lines) calculated using perturbation 
theory. The next-order boundaries for n = 1 are also shown 
(blue dotted lines). 

adding a gain-loss profile that is modulated in space, one 
obtains Eq. Q. Actually, it is not necessary to include 
gain: if there is only modulated loss, one obtains Eq. ([^ 


on top of a background of decay, and the VT -symmetric 
phase corresponds to slower decay [2] . Alternatively, one 
can implement the model with a trapped atom by modu¬ 
lating the spontaneous emission rate (by modulating the 
optical pumping) and post-selecting on trials without de¬ 
cay events [HHIQ]. 

In Sec. [nj we review the perturbation theory for the 
Hermitian Rabi model. In Sec. |III[ we examine the 
non-Hermitian model, considering the single-photon res¬ 
onance, then the multi-photon resonances, and finally the 
limits of low and high frequency. The Appendixes pro¬ 
vide details of the perturbation theory. 


II. REVIEW OF HERMITIAN RABI MODEL 

First, we briefly review the Hermitian Rabi model, 

H = ^az + 2Xax cos wt, (2) 

where the drive strength A is real. We work in the tr^ 
basis: |t), |i). This model was solved perturbatively in 
A in Ref. [20] using the perturbation theory of Salwen 

ESI- 

We move to the Floquet picture (201 EH- We define the 
basis \an) for the Floquet Hamiltonian, where a =t)4- 
and n is any integer. The Floquet Hamiltonian Up is an 
infinite-dimensional matrix given by 12011111 


{an\HF\j3m) = (a| + dm,n-l)\P)- 

For example, the block of Hp corresponding to |t —1), |i —1), It 0), |t 0), |t 1), It 1), It 2), |t 2) is 
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(3) 


(4) 


We seek the eigenvalues e of Hp, since they are the 
quasienergies of H. We are mainly interested in when lo is 
such that two diagonal elements oi Hp are nearly degen¬ 
erate, since this corresponds to a resonance. The main 
(single-photon) resonance occurs when a; « Wq. The 
multi-photon resonances occur when Wq « Sw, Sw, 7w,..., 
i.e., when uto is an odd multiple of ui. 

When two diagonal elements oi Hp are nearly degener¬ 
ate, we use Salwen’s perturbation theory [25] to calculate 
the effective 2x2 Hamiltonian for this subspace. The ef¬ 


fective Hamiltonian captures the effect of all other states 
on the two degenerate states. Once we have this effective 
Hamiltonian, we diagonalize it to obtain e. In Appendix 
[X] we briefly review Salwen’s perturbation theory, since 
it may be unfamiliar to many readers. 
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A. Single-photon resonance: to Ki ljo 


B. Multi-photon resonance: ojo ~ {‘2n+ l)uj 


Suppose that uj ^ ujo, such that the unperturbed states 
It 0) and It 1) are almost degenerate. Using Salwen’s per¬ 
turbation theory (see Appendix]^, we find the effective 
2x2 Hamiltonian H'p for these two states as a pertur¬ 
bation series of A. One must be careful when doing the 
perturbation theory: {uj — uio) and {e — uJo/2) are assumed 
to be on the order of A. 

To lowest order, the effective Hamiltonian is [20] 

= (5) 


Now suppose that ujo ~ (2n -|- l)w where n is a pos¬ 
itive integer. Then the unperturbed states, |t 0) and 
14, {2n + 1)), are almost degenerate. From Salwen’s per¬ 
turbation theory, we obtain the effective 2x2 Hamilto¬ 
nian for these two states: 




f+S u \ 

u + — 5 y’ 


( 12 ) 


where 5 is the level shift, and u is the effective coupling 
between the states. Note that there is no direct coupling 
between |t 0) and |4. (2n -I- 1)) in Eq. Q, so rt arises from 
coupling to intermediate states. The eigenvalues are 


which has eigenvalues. 


' i 


e = 


2 ’ 

= (w - ujof + 4A^ 


tween |4, 1) and |t 2), which are the “counter-rotating” 
terms. 

To next order, the effective Hamiltonian is EDI 



(2n -|- l)w ± H 

2 

(13) 

(6) 

= [(2n -1- l)a; — Wo — 26]“^ + 4u^, 

(14) 

(7) 

One finds [DO] [26] 


e fre- 
r ap- 
since 
d be- 

^ 2n -1- 1 A^ 

2n(n -I- 1) w ’ 

(-1)”A2"+i 

7 / — -- 

22"(n!)2a;2"’ 

(15) 

(16) 


where Eq. (151 is accurate to O(A^), and Eq. (16) is ac¬ 
curate to 0(A^"’+^). Equation shows that the larger 
n is, the weaker the transition, as expected. 






( 8 ) 


III. NON-HERMITIAN RABI MODEL 


where there is now a level shift due to coupling to |4, — 1) 
and It 2). The eigenvalues are now 


w ± H 


e = 


= (w - -I- 4A^ - 


(9) 

( 10 ) 


where H is the effective Rabi frequency. We have trun¬ 
cated at O(A^), because it is only accurate to that 
order. [Improving the accuracy would require calculat¬ 
ing the level shift beyond O(A^), which has been done 
[211126], but we stop here for simplicity.] 


Erom Eq. (10), one finds a shift in the resonance 
frequency Wres, defined as the frequency at which the 
Rabi oscillations have maximum amplitude. Since the 
amplitude of Rabi oscillations is proportional to 1/H^ 
[21j . the oscillations achieve maximum amplitude when 
dfl'^/duj = 0 : 


^res — 


A2 

LOo 


( 11 ) 


Now we adapt the results of Sec. [IT] to the non- 
Hermitian case by replacing A —>■ iA throughout. Al¬ 
though Hp in Eq. 0 becomes non-Hermitian, it is VT- 
symmetric using V = az®I, so we expect the eigenvalues 



This is the famous Bloch-Siegert shift, which is due to 
the counter-rotating terms [DDIED] . 


FIG. 3. (Color online) Imaginary part of Floquet quasienergy 
e for A = 0.2ujo, showing windows of "PT breaking. 
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e to be real for small A and complex for large A. (Al¬ 
though Hp has an infinite number of eigenvalues, they 
all become complex simultaneously, so we only need to 
track one of them.) 

Figures and show that there are windows of VT 
breaking whenever lo hits a resonance. The intuitive rea¬ 
son for this is that when the drive frequency is close to a 
resonance, the imaginary drive has more effect and thus 
breaks VT symmetry. On the other hand, when the drive 
is not close to a resonance, it does not have much effect 
and VT symmetry is maintained. Figure plots the 
imaginary part of e as a function of w. 

We seek to derive the boundary between the VT- 
symmetric and T^T-broken phases. It is clear from 
Eqs. and ( [T^ that e becomes complex when VV 

or rV becomes negative. Thus, the boundary between the 
T^T-symmetric and T^T-broken phases is given by when 
or is zero. 


A. Single-photon resonance: to fa uuo 


We first consider the case of w « Wo, replacing A —>■ zA 
in Sec. II A| We still call this the “single-photon” reso¬ 
nance even though the imaginary drive does not involve 
photons. To lowest order, the eigenvalues are 


' i 


= (w - Wo)^ - 4A^. 

The system is in the T^T-broken phase < 0) when 


(17) 

(18) 


a>F^. (w) 

Thus, the boundary is linear under the rotating-wave ap¬ 
proximation. This was derived previously in Ref. m- 
To next order. 


win , , 

e=^—: (20) 

= (w - Wo)2 - 4A2 i , (21) 

COq 

So the system is in the T^T-broken phase when 

+ ( 22 ) 


4a;„ 


where we have truncated at O(A^), which is the accu¬ 
racy at this level. Thus, we have obtained the quadratic 
correction to the boundary due to the counter-rotating 
terms. This agrees well with numerics for small A (Fig.[^. 

In the non-Hermitian case, there is still a Bloch-Siegert 
shift. When e is complex, the evolution under H leads 
to exponential growth, and the growth rate is maximum 
when reaches its minimum, which happens at 

A2 

a^res “ ■ (23) 

Uq 


At this point, the imaginary part of e reaches its maxi¬ 
mum. Thus, the Bloch-Siegert shift corresponds to the 
point at which the solutions exhibit the fastest exponen¬ 
tial growth. Note that due to the T, the Bloch-Siegert 
shift for the non-Hermitian case is in the opposite direc¬ 
tion of the Hermitian case. 


B. Multi-photon resonance: uio ~ {2n+ l)uj 


At the multi-photon resonances, Wo « (2n -|- l)a;, 
there are windows of VT breaking in the phase diagram 
(Fig. §. Here we obtain the windows to lowest order. 
Adapting the results of Sec. H B to the non-Hermitian 
case: 

(2n -I- l)a; ± H 

= [(2n -|- l)a; — Wo — 2(5]^ -I- 4u^, 

2n -I- 1 A^ 


5 = - 


2n{n -I- 1) w ’ 


(24) 

(25) 

(26) 


22"(n!)^ 


(27) 


We see that at lowest order, does not depend on u. 
However, when we omit u from Eq. (251, cannot be¬ 


come negative; its minimum possible value is zero. This 
means that at lowest order, we can only get a line in pa¬ 
rameter space corresponding to when VT symmetry is 
on the verge of breaking. 

The line is given by when (2n -|- I)w — Wo — 2(5 = 0, or 

1/2 


\{n) = 


n(n -I- l)ijJo 
2n -\-l 


COo 


UJ — 


2n -|- I 


(28) 


where we have used the fact that ojo « (2n -I- I)a;. For 
any n, the line has a square-root shape. As n increases, 
the line becomes steeper. 

This result agrees well with numerics (Fig. |^. For 
small A, each window is very narrow and follows the pre¬ 
dicted line. In fact, Fig.j^a) shows that Eq. ( [^ is useful 
for predicting the frequency where maximal VT break¬ 
ing occurs (where Im e is maximum). In analogy with 
the single-photon case, we label this Wres- 

In order to go further and predict an area for each 
window, we would have to go to high-enough order in 
perturbation theory such that u is included in fV. This 
means we would have to obtain 6 to 0(A'^"), which can 
be tedious. In Appendix [B} we do this for n = I. 

As a rough approximation for the windows, we can 
simply use the leading-order value of 6 given in Eq. (261. 


Although this is not a rigorous perturbative expansion, 
it does give reasonable estimates for the window bound¬ 
aries, which are found by solving (2n + l)a; — Wo — 25 = 
±2iu. We use this approach to estimate the width of the 
window for each n. Approximating nl « •\/27m(n/e)"', 
the window width is 

2(eA)2"+i 


width (^) 


Tm(2n + 1 )^ 0 " 


(29) 
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FIG. 4. (Color online) Properties of multi-photon resonances at A = O.lcao for different n, calculated numerically (black 
circles) and via rough perturbation theory (red crosses), (a) Frequency where maximal VT breaking occurs, (b) Width of VT 
breaking window, (c) Maximum imaginary part of Floquet quasienergy, which determines the maximum growth rate. 


This agrees reasonably well with numerics [Fig. |^b)]. 
The deviations are mainly due to omitting the higher- 
order terms in 5. Thus, as n increases, the width de¬ 
creases exponentially. (We assume A < uJole.) 

We can similarly obtain a rough estimate for the max¬ 
imum imaginary part of e for each multi-photon reso¬ 
nance. This quantity is related to the maximum growth 
rate in the T^T-broken phase. We find 


max Im e(n) 


(eA)^"+^ 
2TTnu}^'^ ’ 


(30) 


which again agrees reasonably well with numerics 
[Fig. El: c)]. So as n increases, the growth rate decreases 
exponentially. This means that in order to observe VT 
breaking for higher n, one must evolve the system for 
longer time (or longer distance in the case of waveg¬ 
uides). Figure El shows example trajectories in the VT- 
broken phase for the single-photon resonance and the 
three-photon resonance (n = 1): the latter clearly grows 
slower than the former. 



Normalized time 


FIG. 5. (Color online) Example solutions in the 'PT-broken 
phase for the single-photon resonance (left) and the three- 
photon resonance (right) with A = O.lWo. The occupation in 
each state is plotted: \c-f\'^ = \ (t| and = | (4-1 


C. Low frequency: cj Uo 

We now consider the limit of low frequency. It is clear 
from Fig.j^that as w —>■ 0, the VT boundary approaches 
A = ojo/i. This can be understood as follows. For small 
w, the non-Hermitian term in Eq. 0 changes very slowly, 
so the Hamiltonian is basically static. We can replace 
2iXcosu}t with its maximum value, 2IA, since VT sym¬ 
metry is most likely to be broken at that point of the 
cycle. Thus, we have 

iJ = -k 2 jA(J2,. (31) 

It is easy to show that this Hamiltonian is in the VT- 
broken phase when A > lOo/4:, in agreement with numer¬ 
ics in the limit w —)■ 0 (Fig. [2]). However, note that if 


one sets A = Wo/d and takes the limit w —>■ 0, one goes 
through an infinite series of multi-photon transitions m- 


D. High frequency: a; 3> oJo 


In the limit of large frequency, the VT boundary con¬ 
tinues to rise (Fig. [^. It is impractical to keep doing 
perturbation theory in A as in Sec. HI A Instead, we can 
treat ujo as the perturbation, and this allows us to obtain 
the boundary in this limit. Here, we just use conven¬ 
tional perturbation theory - not Salwen’s perturbation 
theory. 

[We note that for the Hermitian Rabi model, one typi- 
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cally treats this limit by transforming to a rotating frame 
and then averaging over a period of uj to obtain an ef¬ 
fectively unmodulated system When applied to 

the non-Hermitian model in Eq. Q, this approach pre¬ 
dicts that V7~ symmetry is never broken, which is wrong 
(but see [28]). It turns out that averaging over a period 
does not work for Eq. 0 because it discards information 
about large amplitude oscillations. Thus, we must use 
the alternative approach described below.] 

It is clear from Eq. Q that there are two uncoupled 
groups of states. The hrst group includes |t0), 14-1), 
It 2), It 3). The second group includes |t 0), |t 1), it 2), 
It 3). For the following discussion, we only consider the 
first group. (The second group gives the same results.) 
Eor the first group, we can rewrite Eq. 0 as 

= + (32) 

Hq = -I- I| -I- |m -I- 1 )(to|) -|- maj|TO)(m|], 

m 

(33) 

^1 = ( 34 ) 

m 

where we have replaced A —> iX. Let the eigenvalues and 
eigenstates of Hp be denoted e„ and \4>n)-i where n is 
any integer. (Here, n is unrelated to the n of previous 
sections!) We treat Hi as the perturbation. 

Hq is a well-known Hamiltonian: it is a tight-binding 
model with a linear potential and leads to the Wannier- 
Stark ladder [Mila. The eigenvalues of Hq are 

f-T = nuj, (35) 

and the eigenstates can be written in terms of Bessel 
functions of the first kind, 

= E (~) I"*)- (^®) 

m. ^ ^ 


Note that is always real. 

Now we calculate the first-order shift due to iJi. Eor 
simplicity, we only calculate the shift for the n = 0 eigen¬ 
value: 



which is always real and positive. We have used Bessel- 
function identities in the last two lines. One can show 
that = —Cl. 

However, there is a problem. It is clear that doing per¬ 
turbation theory in this way will never be able to exhibit 
a V7~ breaking transition, since such a transition requires 



cu/o; 

0 


FIG. 6. (Color online) Zoomed-out view of V7~ phase dia¬ 
gram. Also shown is the boundary (red dashed line) calcu¬ 
lated using perturbation theory in cJo. 


two eigenvalues to become degenerate and would thus re¬ 
quire going to infinite order. (Remember that we are not 
using Salwen’s method.) Nonetheless, we can sidestep 
this issue using the following approach. 

We note that the VT transition occurs because as Wo 
increases, eo increases from zero and ei decreases from 
w until the two coalesce and become complex. (Due to 
periodicity of Eloquet eigenvalues, we only have to con¬ 
sider Eq and El.) We can approximate when this degen¬ 
eracy occurs using the fact that eJ°^ — Eg°^ = uj and that 
Eq^^ = — Ej^^ Thus, we expect the degeneracy to occur 
when 2 eq^^ = uj, or 


UJ 

UJo 



(41) 


One has to solve this transcendental equation numeri¬ 
cally. The solution is plotted in Fig. [^ which shows rea¬ 
sonable agreement with numerics. The deviations are 
probably due to the rough way we predicted the degen¬ 
eracy. 


IV. CONCLUSION 

We have derived the VT phase diagram of the 
non-Hermitian Rabi model in different limits. A future 
direction is to extend the results to a many-body setting, 
where the imaginary drive couples to many interacting 
spins instead of a single spin; this may lead to exotic 
magnetic phases li ii mi ED. Also, while we have con¬ 
sidered a classical drive here, it would be interesting to 
consider what happens with a quantum drive, i.e., the 
cos ujt in Eq. 0 would be replaced by creation and an¬ 
nihilation operators m El- The Hermitian quantum 
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Rabi model was exactly solved recently [32], so perhaps 
the solution can be extended to the non-Hermitian case. 
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Appendix A: Review of Salwen’s perturbation 
theory 

Here, we review Salwen’s perturbation theory (Sec. 3B 
of Ref. [21]). Suppose we want to find the eigenvalues 
of a Hamiltonian H, but there are two nearly degenerate 
levels, I a) and \b). Salwen’s method finds the effective 
Hamiltonian H' for the subspace of |a) and |6). This 
effective Hamiltonian is a 2 x 2 matrix that captures the 
effect of all other levels on this subspace. H' is written 
as a perturbation series of the off-diagonal elements. 

Let \ip) be an eigenvector of H corresponding to eigen¬ 
value e: 

H\i:)=e\^). (Al) 

IV') can be expanded as 

IV’) = \a){a\^) + |&)(&|V’) + l”)(’^IV’). (A2) 

n 

where the index n runs over states other than |a) and \b). 
By plugging Eq. (A2) into Eq. (|A1|), we get 


where 


(a|iJ|a)(a|V') + {a\H\b){b\ilj) + ^(a|iJ|n)(n|V') 

n 

= e(alV’). (A3) 

{b\H\a){a\i)) -f (6|i7|6)(6|V') + ^(6|iJ|n)(n|V’) 

n 

= e(&IV’). (A4) 

{n\H\a){a\ip) + (n|i7|5)(&|V’) + ^(n|i7|n')(n'|V') 

n' 

= e(^^IV'),(A5) 


We can solve Eq. (A5| for in terms of {a\'tl)) and 

(6|V’): 


i|V’) = 


A„(e)(a|Vj) -I- Bn{e){b\'ip) 
e — (n|i7|n) ' 


(A6) 


where A„(e) and Bn{e) depend on e and the matrix el¬ 
ements of H. If H is infinite-dimensional (as it is for a 
Eloquet Hamiltonian), then we write A„(e) and Bn{e) as 
power series in the off-diagonal elements of H, keeping 
enough terms to obtain the desired accuracy. 


By plugging Eq. (A6| into Eqs. (A3) and (A4), we 
obtain the effective 2x2 Hamiltonian for |a) and \b): 

_ f {a\H\a}+ Vaa{e) Vab{e) ^ 

' Vba{e) {b\H\b) + Vbb{e)h^^^^ 


Vaaie) = 

n 

Vbbie) = 


{a\H\n)An{e) 
e — (n|iL|n) 

{b\H\n)Bnie) 
e — (nlHln) 




(n\H\n) 


= Vblie) 


(A8) 

(A9) 

(AlO) 

(All) 


In pra ctic e, the procedure is as follows. We write 
out Eq. (A5), then solve it for (n|V’), thereby obtaining 
A„(e) and H„(e) as power series of the off-diagonal ele¬ 
ments to desired order. Then we find H' by calculating 
Vaa, Vbb, Vab- We diagonalize H' to find its eigenvalue e. 
However, if H' itself depends on e, we solve the resulting 
implicit equation for e. 


Appendix B: Boundaries for three-photon resonance 

Here, we calculate the boundaries for the three-photon 
resonance (n = 1) beyond just a single line. This requires 
going to high enough order in perturbation theory. The 
effective Hamiltonian when ojo ~ 3w is: 


H'p = 


-L A 

2 w Oq 


— ^ -f 3w -f 5b 


(Bl) 


where the level shifts of the two states {5a, 5b) can be 
different. The eigenvalues are 


e = 


3w -f l5a -f ^6 


2 ± ^\/4m2 -I- {5a -5b- 3w -I-Wo)2. 


(B2) 

For convenience, we define a small parameter a to do 
perturbation theory in. We rewrite A, e, w in terms of a: 


A = A'a, 

! 2 

e-- = ea , 

a; - 


(B3) 

(B4) 

(B5) 


where A', e', A are on the order of unity. Using Salwen’s 
perturbation theory, we find: 


9A2 2 9A'U9A'2-f bAwo - lOe'wo) 4 

J„ = -—o- 3 - £ ;| .— ° • 

, 9A2 2 9A'2(9A'2-24Aa;o-f lOe'wo) 4 

5b = -— a ^ - 

dWo 


32a;3 


9iA'3 


(B7) 
(B8) 


where 5a,5b are accurate to O(a^). 


We plug Eqs. (B3|“(B8) into Eq. (B2) and solve for 


d. We want to find when d becomes complex, since that 




















corresponds to the VT transition 




0 = (3A'2 + 2/^ujof 

9A'^(5A'^ + 4Aa;o)(llA'^ + 2Awo) 2 

+--a , (B9) 


which can be solved for X' in terms of A. This agrees well 
with numerics for small A (see blue dotted lines in Fig.[2 1. 
Note that if we omitted the second term in Eq. (B9|, 
we would recover Eq. ( [^ , since that is the lowest-order 
result. 


[1] C. M. Bender and S. Boettcher, Phys. Rev. Lett. 80, 5243 
(1998) 

[2] A. Guo, G. J. Salamo, D. Duchesne, R. Morandotti, 
M. Volatier-Ravat, V. Aimez, G. A. Siviloglou, and D. N. 
Ghristodoulides, Phys. Rev. Lett. 103, 093902 (2009) 

[3] E.-M. Graefe, H. J. Korsch, and A. E. Niederle, Phys. 
Rev. A 82, 013629 (2010) 

[4] G. E. Riiter, K. G. Markris, R. El-Ganainy, D. N. 
Gristodoulides, M. Segev, and D. Kip, Nat. Phys. 6, 
192 (2010). 

[5] R. Uzdin, A. Mailybaev, and N. Moiseyev, J. Phys. A 
44, 435302 (2011) 

[6] W. D. Heiss, J. Phys. A 45, 444016 (2012) 

[7] G. Thompson, Y. N. Joglekar, and G. Vemuri, Phys. 
Rev. A 86, 043822 (2012) 

[8] T. E. Lee and C.-K. Ghan, Phys. Rev. X 4, 041001 (2014) 

[9] T. E. Lee, F. Reiter, and N. Moiseyev, Phys. Rev. Lett. 
113, 250401 (2014) 

[10] T. E. Lee, U. Alvarez-Rodriguez, X.-H. Cheng, 
L. Lamata, and E. Solano, Phys. Rev. A 92, 032129 


[11] G. Li and Z. Song, Phys. Rev. A 91, 062104 (2015) 

[12] L.-N. Wu, G.-R. Jin, and L. You, Phys. Rev. A 92, 
033826 (2015) 

[13] M. Znojil, Ann. Phys. 361, 226 (2015) 

[14] J. M. Zeuner, M. G. Rechtsman, Y. Plotnik, Y. Lumer, 
S. Nolte, M. S. Rudner, M. Segev, and A. Szameit, Phys. 
Rev. Lett. 115, 040402 (2015) 

[15] K. V. Kepesidis, T. J. Milburn, K. G. Makris, S. Rotter, 
and P. Rabl, arXiv:1508.00594. 

[16] N. Bender, H. Li, F. M. Ellis, and T. Kottos, 
arXiv:1507.01136. 


[17] X. Luo, J. Huang, H. Zhong, X. Qin, Q. Xie, Y. S. 
Kivshar, and G. Lee, Phys. Rev. Lett. 110, 243902 
(2013) 

[18] Y. N. Joglekar, R. Marathe, P. Durganandini, and R. K. 
Pathak, Phys. Rev. A 90, 040101 (2014) 

[19] J. Gong and Q.-h. Wang, Phys. Rev. A 91, 042135 (2015) 

[20] J. H. Shirley, Phys. Rev. 138, B979 (1965) 

[21] T. Dittrich, P. Hanggi, G.-L. Ingold, B. Kramer, 
G. Schon, and W. Zwerger, Quantum Transport and Dis¬ 
sipation (Wiley-VCH, Weinheim, 1998). 

[22] F. Bloch and A. Siegert, Phys. Rev. 57, 522 (1940) 

[23] I. L. Garanovich, S. Longhi, A. A. Sukhorukov, and Y. S. 
Kivshar, Phys. Rep. 518, 1 (2012) 

[24] K. Shandarova, C. E. Riiter, D. Kip, K. G. Makris, D. N. 
Ghristodoulides, O. Peleg, and M. Segev, Phys. Rev. 
Lett. 102, 123905 (2009) 

]25] H. Salwen, Phys. Rev. 99, 1274 (1955) 

[26] P. K. Aravind and J. O. Hirschfelder, J. Phys. Chem. 88, 
4788 (1984) 

[27] S. Vemuri and Y. Joglekar, (in preparation). 

[28] The averaging method does work when the modulation 
is in the Hermitian term instead of the non-Hermitian 
term m- 

]29] M. Gliick, A. R. Kolovsky, and H. J. Korsch, Phys. Rep. 
366, 103 (2002) 

[30] T. Niemczyk, F. Deppe, H. Huebl, E. P. Menzel, 
F. Hocke, M. J. Schwarz, J. J. Garcia-Ripoll, D. Zueco, 
T. Hummer, E. Solano, A. Marx, and R. Gross, Nat. 
Phys. 6, 772 (2010). 

[31] P. Forn-Dlaz, J. Lisenfeld, D. Marcos, J. J. Garcia-Ripoll, 
E. Solano, G. J. P. M. Harmans, and J. E. Mooij, Phys. 
Rev. Lett. 105, 237001 (2010) 

[32] D. Braak, Phys. Rev. Lett. 107, 100401 (2011) 













